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ABSTRACT 

Context. Although photometric variations of chemically peculiar (CP) stars are frequently used to determine their rotational periods, 
the detailed mechanism of their light variability remains poorly understood. 

Aims. We simulate the light variability of the star HR 7224 using the observed surface distribution of silicon and iron. 
Methods. We used the TLUSTY model atmospheres calculated for the appropriate silicon and iron abundances to obtain the emergent 
flux and to predict the rotationally modulated light curve of the star. We also obtained additional photometric measurements and 
employed our own regression procedure to derive a more precise estimate of the light elements. 

Results. We show that the light variation of the star can be explained as a result of i) the uneven surface distribution of the elements, 
ii) the flux redistribution from the ultraviolet to the visible part of the spectrum, and iii) rotation of the star. We show that the silicon 
bound-free transitions and iron bound-bound transitions provide the main contribution to the flux redistribution, although an additional 
source of opacity is needed. We confirm that numerous iron lines significantly contribute to the well-known depression at 5200 A and 
discuss the connection between iron abundance and the value of peculiarity index a. 

Conclusions. The uneven surface distribution of silicon and iron is able to explain most of the rotationally modulated light variation 
in the star HR 7224. 

Key words, stars: chemically peculiar - stars: early type - stars: variables - stars: individual HR 7224 



1. Introduction 

Chemically peculiar (CP) stars represent a large class of upper 
main sequence stars where the processes of radiative diffusion 
and gravitational settling in their atmospheres give rise to pro- 
nounced deviations in the chemical composition of these stars 
from the solar value (Vauclair 2003, Michaud 2005). The CP 
stars are natural laboratories for testing modern model atmo- 
spheres thanks to the unusual chemistry with rather strong under- 
or overabundance of some elements. The application of ad- 
vanced modelling techniques, such as model atmospheres with 
magnetic fields (e.g., Kochukhov et al. 2005, Khan & Shulyak 
2006), radiative diffusion codes, and Doppler imaging tech- 
niques (Khokhlova et al. 2000), provides us with detailed infor- 
mation about the surface structure of these stars. Despite these 
fascinating advances in their study during recent decades, the 
light variability of CP stars is still poorly understood. 

Some CP stars show inhomogeneous surface distribution of 
chemical elements on their surface (e.g., Khokhlova et al. 2000) 
as determined from rotationally modulated spectral line variabil- 
ity (see, e.g., Lehmann et al. 2006). The uneven surface distri- 
bution of various elements, together with rotation, has been pre- 
sumed to be the origin of these stars' light variability. However, 
the details of this mechanism have not been determined. Line 
blanketing by multiple lines of overabundant elements (mainly 
iron) and the flux redistribution induced by these lines has been 
proposed as one of the causes of the light variability (Molnar 
1973). Other mechanisms proposed include the influence of 
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bound-free transitions (Peterson 1970; Lanz et al. 1996), sur- 
face temperature differences or variable temperature gradients 
(Kodaira 1967; Weiss et al. 1976), and the presence of magnetic 
fields (e.g., Kochukhov et al. 2005). Finally, circumstellar mat- 
ter, if present, may also influence the light curves (Landstreet & 
Borra 1978; Nakajima 1985; Smith & Groote 2001; Townsend 
et al. 2005). 



One of the first attemps to simulate the light variability of 
CP stars was done by Krivosheina et al. (1980), who reproduced 
the light curve of the CP star CU Vir. However, detailed mod- 
elling of CP star light variability had to await precise model at- 
mospheres (e.g., Lanz & Hubeny 2007), improved opacity data 
(Seaton et al. 1992), and much faster computers. Krticka et al. 
(2007) took advantage of these tools and used the surface maps 
of Khokhlova et al. (2000) to simulate the light curve of HD 
37776 successfully. They demonstrate that the inhomogeneous 
surface distribution of silicon and helium, along with the bound- 
free transitions of these elements, accounted for most of the light 
variability in this star. 



Much work remains to be done to understand the light vari- 
ability in CP stars. In particular, the role of iron, which is 
found to be significantly underabundant in the atmosphere of HD 
37776 (Khokhlova et al. 2000), remains to be clarified. For this 
purpose we selected the silicon star HR 7224, whose silicon and 
iron surface distributions were derived via Doppler imaging by 
Lehmann et al. (2007). 
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2. CP star HR 7224 

HR 7224 (HD 177410, HIP 93187, EE Dra) is a rapidly rotat- 
ing CP star classified as AOp Si by Cowley et al. (1969). Molnar 
(1972) classified HR 7224 as a B9.5 IIIp Si star, based on low- 
resolution spectroscopy. Lehmann et al. (2006, 2007) found that 
the star has enhanced Si and Fe abundances, whereas He is ex- 
tremely depleted; they classified it as a helium-weak silicon CP 
star. 

Winzer (1974) found photometric variability in HR 7224 
and determined a period of l'?1663 from 48 UBV observa- 
tions done in 1970-72. Adelman (1997) derived a new period 
of l'?123095 based on 616 Stromgren uvby observations taken 
in the 1993-94 and 1994-95 observing seasons with the Four 
College Automated Photometric Telescope (FCAPT) along with 
the older V measurements of Winzer (1974). Hipparcos obtained 
409 Hp, Bj, and Vt observations in 1989-92 (ESA 1997), and 
ESA (1998) determined the photometric period to be l'?123248. 

All of the optical light curves phase together fairly well and 
show a double-wave with two maxima of different height (see 
Fig. 1). HR 7224 was recently classified from the shape of its 
light curve as a prototype of a photometrically simple CP star 
with a double-wave light curve with two unequally prominent 
bright spots centred on the phases (p - 0.0 and 0.5 (Mikulasek 
et al. 2008a). There are only two attempts to measure the mag- 
netic field of HR 7224 available in the literature (Bohlender et 
al. 1993; Lehmaim et al. 2007); both produced negative results. 
Recent magnetic field measurements (Kudryavtsev, private com- 
munication) also gave negative results, however the observations 
are still continuing to cover the whole rotational cycle. 

Adelman (2004) reported an unprecedented change in the 
photometric behaviour of HR 7224. Comparing time series ob- 
servations taken with FCAPT before 1996 with those taken in 
2003, he found that the amplitude of variability had increased 
from a typical value of 0.04 mag to 0.21 mag. He also reported 
that the period of variability had changed from l'?123 to 101 
days. Unfortunately, the photometric data from this critical era 
have not been published. 

Adelman's astonishing result, along with the lack of avail- 
able high-resolution spectroscopy of HR 7224, motivated 
Lehmann et al. (2006) to carry out an extensive spectroscopic 
observing program on the star. Their 564 high-resolution spec- 
trograms allow determination of radial velocities via cross cor- 
relation to an accuracy better than 100 m s"^ . They found radial 
velocity variations in this swiftly rotating star with an amplitude 
of 15 km s~^ and a period of if 123248(9), in excellent agreement 
with the previous determination by ESA (1998). They found no 
further periodicities, in particular nothing around the 101-day 
period of Adelman (2004). 

Lehmann et al. (2007) used the same 564 spectrograms to 
derive Doppler images of surface elemental distributions. Their 
map of the silicon and iron abundance on the surface of HR 7224 
is an ideal starting point for the simulation of the star's expected 
light curves. To avoid any difficulties phasing the spectroscopic 
and photometric observations together, we acquired new pho- 
tometric observations to improve the rotational ephemeris of 
HR 7224. 

2.1. New BV photometry of HR 7224 

Our new BV photometry was acquired between March and May 
2008 with the T3 0.4 m automatic photoelectric telescope (APT) 
at Fairborn Observatory. This APT uses a temperature-stabilised 
EMI 9924B photomultiplier tube to measure photon count rates 
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Fig. 1. Photometric observations of HR 7224 plotted as a function of 
the new linear phase in Eq. (3). Data are from Winzer (1974), BV (A); 
Adelman (1997), uvby (□); Hipparcos (ESA 1997), Hp (o), and the T3 
APT (this paper), BV (❖). Solid lines denote the fit according to Eq. (1). 
BtVt and U magnitudes are not plotted here due to their large scatter. 



through Johnson B and V filters. HR 7224 and its comparison 
stars were measured in the following sequence, termed a group 
observation: K-S-C-V-C-V-C-V-C-S-K, where K is the check 
star (HD 172728, V = 5.74, B-V ^ -0.045, AO V), C is the com- 
parison star (HD 172569, V = 6.07, B-V = 0.279, FO V), V is 
HR 7224, and S is a sky measurement. Three V-C and two K-C 
differential magnitudes are formed from each sequence and aver- 
aged together to create group mean differential magnitudes. The 
typical precision of a group mean is 0.003-0.004 mag for this 
telescope. To filter out observations taken in non-photometric 
conditions, group means with a standard deviation greater than 
0.01 mag were discarded. From the 572 group observations col- 
lected by the telescope, 478 B and 462 V group means survived 
the filtering process and were used in this analysis. 

The surviving group means were corrected for differential 
extinction with nightly extinction coefficients, transformed to the 
Johnson system with yearly-mean transformation coefficients, 
and treated as single observations thereafter. Up to five group 
observations were acquired each clear night at intervals of ap- 
proximately two hours. Further information on the operation of 
the APT and the analysis of the data can be found in Henry 
(1995a,b) and Eaton, Henry, & Fekel (2003). 

The photometric data used in this analysis are available 
through SIMBAD and the On-line database of photometric ob- 
servations ofmCP stars (http://astro.physics.muni.cz/mcpod). 

2.2. New linear ephemeris 

To derive the ephemeris, we used all available photometric 
data on HR 7224 including the UBV measurements of Winzer 
(1974), the Hp and BjVj observations from Hipparcos (ESA 
1997), the Stromgren uvby photometry by Adelman (1997), and 
the T3 APT BV observations in this paper. The total of 2013 
individual photometric measurements covers a time span of 38 
years; they were supplemented with the 564 measurements of 
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silicon line equivalent widths (EW) obtained by Lehmann et al. 
(2007) in 2004-05 and kindly provided to us by Dr. H. Lehmann. 

The analysis of these data to determine the ephemeris is de- 
scribed in Mikulasek et al. (2008b). It is assumed that all the 
photometric and spectroscopic variability can be expressed in 
terms a simple regression model. 

2.2.1. Regression model and robust regression procedure 

The formulation of the observed phase variation of the light and 
equivalent widths of silicon lines in terms of a general regres- 
sion model is relatively simple since the light and EW phase 
curves are similar in shape, though shifts may exist in some 
cases. Consequently, we can express them with the following 
relations: 



ycjiti) 



Fi&ij), &ij = {ti-Mo-AMoj)/P, (1) 



where ycj is the predicted value of an observed quantity 
(magnitude, EW) in a colour (pass band) c (in this case 
M, U, V, B, b. Hp, y, V, and EW) of the j-th set of measure- 
ments (Winzer, Hipparcos, Adelman I and II, Lehmann, and this 
paper) at the particular JDhei instant and y^j is the weighted- 
mean value of the quantity. Note that observations in Bt were in- 
cluded into B and Vt into V. Here Ac is an "effective amplitude", 
the parameter robustly describing the measure of the variability 
in the colour c, and F(i?) is a normaUsed (unique) function ex- 
pressing the form of the observed phase curves (for details see 
Mikulasek et al. 2007a; Mikulasek et al. 2007b). It is a periodic 
function of the monotonically growing "phase function" i? intro- 
duced by Mikulasek et al. (2008b) 

The phase function of time §{t) is determined by the 
ephemeris, which is given in the second part of Eq. (1). We as- 
sume that the period P of HR 7224 is constant and define zero 
phase to be the light maximum Mq that is situated nearest to 
the weighted centre of gravity of measurements used for the 
ephemeris determination. AMoj are free parameters correcting 
Mo for some groups of observation. We applied them only when 
justified, particularly for observations of Winzer (1974) (see the 
following discussion) and for the equivalent widths of the Si 
lines. 

Function F(i?) is the simplest normalised periodic function 
that represents the observed photometric and spectroscopic vari- 
ations of HR 7224 in detail. The phase of maximum brightness 
is defined to be 0.0, and the amplitude is defined to be 1 .0. The 
function, being the sum of three orthogonal terms, is described 
by two parameters a\ and a2. The first parameter quantifies the 
symmetric portion of the deviation of the light curve from a sim- 
ple cosine course; the second parameter expresses any asymme- 
try in the light curve: 



F(§, a\, a2) = yjl-a^-al cos(27r?5-) -i- ai cos(47r)9') 
sin(2 n&)-^ sin(4 n §) 



+ a2 



(2) 



2.2.2. New findings 

We used the same robust, iterative, regression procedure de- 
scribed in Sect. 3.3 of Mikulasek et al. (2008b) to determine 
a new linear ephemeris for HR 7224. The ephemeris consists of 
the period P, the JDhei of the light curve maximum Mq, the mean 
shift of Winzer's measurements AMowin, and the mean shift of 
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Fig. 2. 0-C diagram of HR 7224 computed with the new linear 
ephemeris in Eq. (3). The open star represents the observations of 
Winzer, the filled diamond corresponds to the Hipparcos measurements, 
the filled circles are from the two seasons of observations by Adelman, 
and the filled square gives the result of the BV observations in this pa- 
per. The solid line denotes the line of constant period and the dotted line 
is the uncertainty. 



the maximum of the Si equivalent width, AMqew, relative to the 
light-curve maxima: 



P = l'?1232524(6), 
AMowin = O?058(10), 



Mo = 2451582.7164(12), 
AMoEw = 0?0341(32). 



(3) 



Table 1. O-C derived from photometry and spectroscopy of HR 7224. 
s denotes standard deviation in mmags. 



Year 


O-C [d] 


s 


N 


Data 


Source 


1972 


0.0580(70) 


4.4 


48 


UVB 


Winzer (1974) 


1991 


-0.0010(60) 


11 


409 


BtVjHp 


ESA (1997) 


1994 


-0.0005(20) 


3.8 


317 


uvby 


Adelman (1997) 


1995 


0.0005(18) 


3.3 


299 


uvby 


Adelman (1997) 


2006 


0.0341(24) 




564 


EW 


Lehmann et al. 


2008 


0.0000(17) 


4.8 


940 


BV 


this paper 



Our period agrees with those determined by ESA (1998) 
and Lehmann et al. (2006) within their uncertainties; the period 
appears to be constant over the past 17 years (see Fig. 2 and 
Table 1). The nonzero shift AMoew is likely connected with the 
fact that uneven horizontal distribution of silicon is not a unique 
source of the light variabihty (see Sect. 6.1). Winzer's (1974) 
UBV light curves have the same shape as the more recent ones 
but are shifted with respect to the new ephemeris by O?058(10) 
(6cr). Accepting the reality of the shift, we can speculate about 
the possibility of a transient rotational braking that might have 
occurred before the Hipparcos mission. We discussed a sim- 
ilar event in the He-strong star HD 37776 (Mikulasek et al. 
2008b). In the case of HR 7224 however, we have no obser- 
vational evidence of such deceleration between the Winzer and 
Hipparchos epochs. If we assume a constant rate of deceleration 
{P = const.), then we need P ^ 3.4 x 10"^ in order to explain 
the of O-C shift of Winzer's observations. Such a deceleration 
would manifest itself in the lengthening of the period P by 1.8 
second during this time interval of 17 years (1972 - 1989). 
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Fig. 3. Shapes of the normahsed hght curves of HR 7224 before 2000 
(sohd hnes) and after 2000 (dashed hnes). The outer light lines demar- 
cate the uncertainty in the course of the light curves. 



The parameters ai, 02, which describe the shape of the nor- 
malised function F(§, ai, ai) (Eq. (2)) show some minor differ- 
ences before and after the year 2000: 

before 2000 ai = 0.610(10); 02 = 0.084(12); 
after 2000 ai = 0.560(13); ^2 = 0.018(15). (4) 

If real, this means that the secondary maximum of HR 7224 has 
been decreasing in brightness while the two minima have be- 
come nearly identical. It should be noted that the time depen- 
dence of parameters ai, ai has no obvious effect on the star's 
period. 

Changes in the shape of light curves are very rare among 
CP stars: the only star with well-documented light curve varia- 
tions is the rapidly rotating silicon CP star 56 Arietis. The light 
curve variations are accompanied by a steady lengthening of 
the stellar period (Ziznovsky et al. 2000; Adelman et al. 2001). 
On the other hand, in the most rapidly braked CP star known 
(HD 37776), we did not find any secular changes in its hght 
curve (Mikulasek et al. 2008b). However, it is premature to spec- 
ulate further on the nature of long-term light-curve variations of 
HR 7224 without further photometric and spectroscopic moni- 
toring. 

To compare the observed and calculated light-curve shapes 
and phases, we determined the shift between phases computed 
with our ephemeris and the ephemeris used by Lehmann et 
al. (2007) for surface mapping. The observed rotational vari- 
ations of the silicon Une equivalent widths demonstrated that 
Lehmann's phases precede ours by 0.1955(25). We have cor- 
rected for this phase difference in our analysis of the elemental 
surface distributions. 



3. Calculation of the light curve 

3.1. Stellar parameters 

The stellar parameters of HR 7224 adopted in this study, as well 
as the abundances of silicon and iron derived from the abun- 
dance maps (see Fig. 4) and the abundance of helium are taken 
from Lehmann et al. (2006, 2007, see Table 2). The abundances 
are relative to hydrogen, eei = log (A^ei/A^n)- Note that the abun- 
dances used by Lehmann et al. (2006, 2007) are defined slightly 
differently as log(Afei/A'totai)- 



It vanabihty of the Si star HR 7224 
Table 2. Stellar parameters of HR 7224. 



Effective temperature Teff 


14 


500 K 


Surface gravity log g (CGS) 




4.2 


Inclination / 




65° 


Hehum abundance 




= -3.5 


Sihcon abundance 


-3.6 < 


£si :£ -2.4 


Iron abundance 


-4.4 <, 


epe < -3.0 



3.2. Model atmospheres and synthetic spectra 

We used the code TLUSTY (Hubeny 1988; Hubeny & Lanz 
1992, 1995; Lanz & Hubeny 2003) for the model atmosphere 
calculations. We computed plane-parallel model atmospheres in 
LTE, taking into account the atomic data suitable for B type 
stars (Lanz & Hubeny 2007). The atomic data for silicon are 
taken from Mendoza et al. (1995), Butler et al. (1993), and 
Taylor (2008), in preparation; for iron from Kurucz (1994), 
Nahar (1996), Nahar (1997), Bautista & Pradhan (1997), and 
Bautista (1996), and for other elements from Luo & Pradhan 
(1989), Fernley et al. (1999), Tully et al. (1990), Peach et al. 
(1988), Hibbert & Scott (1994), and Nahar & Pradhan (1993). 

We assumed fixed values of the effective temperature and 
surface gravity and adopted models with a generic value Vtmb = 
2kms"'. We assumed the number density of helium relative to 
hydrogen to be log (NhJNh) = -3.5 (Lehmann et al. 2007). The 
abundance of sihcon and iron differed in individual models as 
explained below. We used the solar abundance of other elements 
according to Asplund et al. (2005). 

For the spectrum synthesis (from which we calculate the 
photometric colours), we used the SYNSPEC code. We took into 
account the same transitions as for the model atmosphere cal- 
culations, and we also included the lines of all elements with 
the atomic number Z < 30 not included in the model atmo- 
sphere calculation. We computed angle-dependent intensities for 
20 equidistantly spaced values of fi - cos 6, where is the angle 
between the normal to the surface and the line of sight. 

The model atmospheres and angle-dependent intensities 
I{A, 9, Esi, fipe) were calculated for a two -parametric grid of sili- 
con and iron abundances £si = log (Nsi/Nu) = -4.0, -3.5, -3.0, 
-2.5, and -2.0, and epe = log (A^e/A^n) = -4.5, -4.0, -3.5, and 
-3.0. 

3.3. Phase dependent magnitude 

The radiative flux observed at the distance D from the star with 
radius /?» in a colour c is calculated as 

/e = (^) j /c(e,n)cos0dQ, (5) 

visible 
surface 

where the intensity lc(0, Q.) at each surface point with surface 
spherical coordinates £2 is obtained by means of interpolation 
between the intensities^ IdO, £si. spe) calculated from the grid of 
synthetic spectra as 

^c{A)I{A,0,esi,ep,)dA. (6) 





' Note that the SYNSPEC code outputs the specific intensities per 
unit frequency (not per unit wavelength). 
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Fig. 4. Observed distribution of silicon (left) and iron (right) on the visible surface of HR 7224 at different rotational phases (Lehmann et al. 2007). 
Phases were calculated with the parameters given in Eq. (3). 



The transmissivity function of a given filter c of the 

Stromgren and Aa photometric systems is approximated for sim- 
plicity by a Gauss function that peaks at the central wavelength 
of the colour Ac with dispersion cr^. (see Table 3). The values for 
the uvby (Stromgren) photometric system are taken from Cox 
(2000). Values for the Aa system (colours gi and g2) are taken 
from Kupka et al. (2003). 



Table 3. Central wavelengths and dispersions of the Gauss filter simu- 
lating the transmissivity functions. 



Colour 


u 


V 


b 


Si 


82 


y 


^ [A] 


3500 


4100 


4700 


5000 


5220 


5500 




230 


120 


120 


65 


65 


120 



The observed magnitude difference is 



Arric = -2.5 log ( — 

Tc 



(7) 



where fc is calculated from Eq. (5) and fl'^^ is the reference flux 
obtained under the condition that the mean magnitude over all 
phases is zero. 

4. Influence of abundance anomalies on emergent 
fluxes 

As shown in Fig. 5, the enhanced silicon abundance results in 
the increase of the temperature in the continuum-forming re- 
gion (within Tioss ~ 0.1 - 1) of the model atmosphere. The 
increased temperature is caused by enhanced silicon opacity in 



the model atmosphere. To understand the contribution of bound- 
free and bound-bound (line) processes to the temperature in- 
crease, we calculated model atmospheres with enhanced sili- 
con abundance but neglecting silicon line transitions^ (see also 
Fig. 5). Neglecting silicon line transitions does not significantly 
influence the atmospheric temperature; consequently, enhanced 
opacity due to silicon bound-free transitions is the main cause of 
the temperature increase. 

The silicon bound-free transitions are important mainly in 
the ultraviolet (UV) spectral regions at wavelengths shorter than 
1600 A. In atmospheres with overabundant silicon, the short- 
wavelength part of the spectrum is redistributed to the longer 
wavelengths of the UV spectrum, and also to the visible spec- 
tral regions (see Fig. 6). Consequently, the silicon-rich spots are 
bright in the uvby colours, and are dark in the ultraviolet bands 
withi < 1600 A. 

A similar situation occurs for iron overabundance. Enhanced 
iron abundance leads to an increased temperature in the 
continuum-forming region between T,oss ~ 1 andTi-oss ~ SxlO"-' 
(see Fig. 5). Unlike silicon, enhanced iron abundance warms the 
atmosphere above Tross ^ 5 x 10"^ mainly by line transitions. 
On the other hand, the outermost regions (Tross S 10 '*) of the 
model with spe = -3 are slightly cooler than those with lower 
iron abundance. This effect is also due to iron lines, as the tem- 
perature is higher in the model with neglected line opacity (see 
also Fig. 5). A similar effect was reported by Khan & Shulyak 
(2007). 

The influence of iron overabundance on the spectrum is 
more complicated than the influence of silicon overabundance, 
as shown in Fig. 7. There are several depressions in the UV 



- Note that the resonance structure of silicon bound-free transitions 
was included in this test. 
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Fig. 5. The dependence of temperature on the Rosseland optical depth Tross in atmospheres with various chemical compositions. Left: Influence 
of silicon abundance on the temperature for = -4.5. Right: Influence of iron abundance on the temperature for ggi = -4.0. Crosses denote 
atmospheres with enhanced silicon or iron abundance, but neglecting the opacity due to line transitions of these elements. 
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Fig. 6. The emergent flux from atmospheres with different silicon abundances. The flux was smoothed by a Gaussian filter with a dispersion of 
10 A to show the changes in continuum, which are important for photometric variabihty. The passbands of the uvby photometric system are also 
shown on the graph (gray areas). The fluxes were calculated with TLUSTY for an iron abundance epe = -4.5. 




Fig. 7. The same as Fig. 6 but for iron. The fluxes were calculated with TLUSTY for a silicon abundance ssi = 



-4.0. 
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Fig. 8. The magnitude difference Am^ (see Eq. (8)) between the emer- 
gent fluxes calculated with enhanced abundance of either siUcon or iron 
and fluxes calculated assuming ssi - -4 and gpe = -4.5. The fluxes 
were heavily smoothed by a Gaussian filter with a dispersion of 100 A. 
The depression at 5200 A due to numerous iron lines is clearly apparent. 



spectral region caused by numerous iron lines, but the flux with 
A > 2500 A increases with increasing iron abundance. Generally, 
the iron hne transitions redistribute the emergent UV radiation 
with A < 1700 A primarily to the long wavelength part of the 
spectrum with A > 2500 A. Consequently, iron-rich regions are 
bright in the uvby colours, whereas they are dark in the ultravio- 
let bands with A < 1700 A. 

These flux changes can be detected as a change in the appar- 
ent magnitude of the star (see Fig. 8). Here we plot the relative 
magnitude difference defined as 



Am^ = -2.5 log 



'//^(esi.epe) 



(8) 



where IPf is the reference flux for ssi - -4.0, spe = -4.5. For 
both silicon and iron overabundance models, the absolute value 
of the relative magnitude difference decreases with increasing 
wavelength. For the silicon model, the decrease is nearly fea- 
tureless, whereas the iron model exhibits several depressions in 
Am^ due to the cumulation of iron lines. The most prominent 
feature at about 5200 A contributes significantly to the 5200 A 
depression frequently observed in the spectra of CP stars. Khan 
& Shulyak (2007) arrived at the same conclusion. 

5. Predicted light variations 

Predicted light curves are calculated from the surface abundance 
maps derived by Lehmann et al. (2007) and from the emergent 
fluxes computed with the SYNSPEC code, applying Eq. (7) for 
individual rotational phases. Since Lehmann et al. (2007) pro- 
vide surface maps with different regularisation parameters A, we 
selected a map with A = 0.0001, which is the most detailed one, 
for our initial calculations. 

To study the influence of silicon and iron separately, we first 
calculated the light variations due to silicon only, assuming a 
constant iron abundance of epe - -4.5. The observed light max- 
imum occurs at the same phase at which the silicon Unes have 
their maximum strength {(/> = 0). As sihcon rich regions are 
bright in visible bands, our predicted light maximum should also 
occur at this phase. Indeed, there is a good agreement of both 
the times of maxima and minima of our predicted light curve 



0} 
•D 




0.2 0.4 0.6 0.8 

phase (t> 

Fig. 9. Light variations of HR 7224 calculated from the silicon abun- 
dance maps only (dashed line) and from the iron abundance maps only 
(solidUne). For both cases, A = 0.0001. Observed light variations (open 
circles) are taken from Adelman (1997). 



and the observed light curve from Adelman (1997), though the 
amplitudes are different (Fig. 9). 

A similar test was performed using only the iron abundance 
map. The silicon abundance was assumed to be esi = -4.0. 
The iron lines are also observed to have their maximum strength 
close to phase = 0, i.e., during observed light maximum. As 
the iron-rich regions are bright in the uvby colours (see Fig. 7), 
the predicted hght curve due to iron abundance variations alone 
also has a hght maximum at phase = 0, in agreement with the 
observations. The predicted light curves due to iron only, dis- 
played in Fig. 9, also have lower amplitudes than the observed 
ones. 

Including the surface distribution of both sihcon and iron in 
the calculation of the light curves (Fig. 10), we obtain very good 
agreement between the observed and predicted light curves in 
the V, b and y wavelength bands. In the u band, however, the 
observed maximum is higher and the first minimum is deeper 
than predicted. A similar situation is seen in the colour indexes: 
(v - b) and (b - y) agree well with the observations while the 
(m - b) curve is not a good fit, due to the differences in the u 
band fit (Fig. 11). The observed Balmer discontinuity index ci 
exhibits large scatter due to the combination of the noisy u, v 
and b light curves, but its variation is nonetheless in accordance 
with the observations of Lehmann et al. (2007) on a possible 
variability of H/3 line. Similarly, the mi index, which reveals the 
effect of the metalhc bound-bound absorption in the violet when 
compared with the b and y regions, shows a phase-dependence. 
The differences between the predicted curves of c\ and mi and 
the observed ones clearly point to the existence of an additional, 
unknown mechanism working in the violet band that still has to 
be investigated (Fig. 10, and see also Sect. 6.1). 

Thus, the inhomogeneous surface distribution of silicon and 
iron revealed in spectroscopic observations results in the appear- 



8 



J. Krticka et al.: The nature of the hght variabihty of the Si star HR 7224 




-0.02 



0.2 0.4 0.6 0.8 

phase 4> 

Fig. 10. Predicted light variations of HR 7224 computed taking into 
account both the silicon and iron surface abundance distributions de- 
rived by Lehmann et al. (2007) for different regularisation parameters 
A. Observed light variations (circles) are taken from Adelman (1997). 



ance of spots on the stellar surface causing the photometric vari- 
ability (see Fig. 12). 

6. Discussion 

6.1. Detailed comparison of observed and predicted liglit 
curves 

The predicted and observed light curves differ slightly in their 
shapes (Fig. 10) and amplitudes (Fig. 13). The small differences 
in the observed vs. calculated amplitudes, leads to the differences 
in the observed vs. calculated colour indexes, especially for ci 
and nil (see Fig. 11). The difference between the observed and 
predicted light curves in Fig. 10 could be explained by the pres- 
ence of an additional photometric spot(s) on HR 7224 caused 
by overabundance of other element(s). We suspect chromium or 
magnesium. Khan & Shulyak (2007) suggested chromium as an 
element capable of significant influence over the emergent flux 
from a star. On the other hand, Lehmann et al. (2006) found 
variability in the magnesium lines of HR 7224 with maximum 
strength occurring roughly at the same phase as the maximum 
difference between theory and observation in Fig. 10 (slightly 
before the maximum strength of silicon Unes). This, along with 
the fact that the calculated effect of the magnesium abundance 
on the emergent flux has a significantly larger amplitude in the 
u band rather than in v, b and y, argues for magnesium as the 
cause of the residuals between observed and calculated light 
curves when considering only silicon and iron. Note however, 
that the maximum magnesium abundance derived by Lehmann 
et al. (2006) without surface mapping is too small to affect the 
spectral energy distribution significantly. The same is also true 
for helium and oxygen. 

Lehmann et al. (2007) derived three abundance maps calcu- 
lated with a different regularisation parameter. A, correspond- 
ing to a different assumed minimum size of the surface inho- 
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Fig. 11. Predicted variations of colour indexes (soUd line) calculated 
from the surface abundance maps derived with A = 0.0001 compared 
with the observations. 



mogeneities. The shapes of the light curves calculated using the 
abundance maps obtained with a different value of A are similar 
(see Fig. 10). A detailed comparison of these curves shows that 
better agreement between theory and observation is obtained for 
the most complex surface map, i.e., for A = 0.0001. However, 
the difference between observed and predicted light curves is 
significantly larger than the difference among the individual pre- 
dicted curves calculated with different values of A. 



6.2. Modei assumptions 

There are several effects that can influence the predicted light 
curves (see Krticka et al. 2007, for a more detailed discussion). 
For example, NLTE effects may influence the continuum flux 
distribution, but since HR 7224 surface maps were derived using 
LTE models, we confine ourselves to LTE models only. Strong 
surface magnetic field also influences the emergent spectral en- 
ergy distribution (e.g.. Khan & Shulyak 2006). However, as all 
available measurements of surface magnetic field are negative 
and the corresponding upper limit is too low to signicantly in- 
fluence the emergent flux, we neglect the influence of magnetic 
field. 

HR 7224 is classified as a helium-weak chemically peculiar 
star, so a question could arise concerning the effect of the uneven 
surface distribution of helium on the light curve. In Krticka et al. 
(2007) we showed that helium influences the spectral energy dis- 
tribution only in the case when it is overabundant (with respect 
to solar), so we neglect any possible inhomogeneous helium sur- 
face distribution in our model. We also tested how accumulation 
of He in the sub-photospheric layers affects emergent flux and 
draw a conclusion that for Tross > 1 this has no effect on the light 
variabihty. 
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Fig. 12. The emergent intensity (in the y band, yu = 1) from individual 
surface elements of HR 7224 at various rotational phases. 
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Fig. 13. The effective amplitudes of the predicted and observed light 
curves plotted versus the effective wavelengths of the various bands. 
The observed light-curve amplitudes are plotted with filled circles (•). 
The predicted amplitudes calculated using both silicon and iron maps 
are filled diamonds (♦), amplitudes calculated using iron only are open 
squares, (□), silicon only are triangles (A), and finally the residuals be- 
tween observation and theory as open circles (o). Observed amplitudes 
in BtVt and U are not plotted here due to their large uncertainty. 



6.3. UV variations 

The redistribution of the flux from short to longer wavelengths 
is one of the consequences of the proposed mechanism for light 
variability in HR 7224. From our calculated model fluxes in 
Figs. 6, 7, we predict that the star should have a flux minimum 
in the short-wavelength A < 1600 A part of the spectrum during 
the visible light maximum and vice versa. A similar antiphase 




-0.06 -0.05 -0.04 -0.03 -0.02 -0.01 0.01 
S {b-y) [mag] 

Fig. 14. The dependence of the Sa index on S(b - y) for individual mod- 
els from the model grid. The Sa values for models with the same iron 
abundance are connected with solid lines, and Sa values for the models 
with the same silicon composition are connected with dashed lines. 




Fig. 15. The light curve in Sa index 

behaviour of the short-wavelength UV and optical light curves 
has been reported for other CP stars (e.g., Sokolov 2006). 

Unfortunately, there is no UV light curve of HR 7224 avail- 
able. However, the star was observed in UV by the TDl space- 
craft (Jamar 1977, 1978). We note that our fluxes calculated with 
overabundant iron in Fig. 7 are in qualitative agreement with 
those observed by TDl at dififerent phases. In accordance with 
observations, the maximum UV light amplitude occurs at about 
1400 A, and the variations at this wavelength are anticorrelated 
with the variations in the visible. The variations at 1600-1900 A 
are very small, and the variations at 2100 - 2300 A are inversely 
correlated with that at 1400 A. The absence of observed light 
variations at 2740 A may be connected with strong iron lines in 
this spectral region. Finally, the observed amplitude of the light 
variations at 1400 A (which is 0.24 mag, Jamar 1978) is in very 
good agreement with those derived from our models (0.21 mag 
for a Gaussian filter with dispersion 50 A). 

6.4. Variations of a 



The peculiarity index a introduced by Maitzen (1976) as 
a^gi- igi +y)/2 



(9) 



is based on the existence of flux depressions (Kodaira 1969) 
correlated with the chemical peculiarity of a given star To test 
whether the peculiarity demonstrated in the a-index is indeed 
caused by numerous iron lines (Khan & Shulyak 2007), we cal- 
culated the peculiarity index a for individual model atmospheres 
used in our previous calculations. 

The dependence of a on the iron and silicon abundance can 
be seen in Fig. 14, where we plot the dependence of 6a on 6(b - 
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Fig. 16. The calculated variations of H6 Une depth with phase 

y). These indices are defined relative to their values for the model 
with esi = -4.0 and epe = -4.5, 

6a =a - aesj=_4.o,6p,=-4.5, (10) 
6{b - y) ={b -y)-(b- y).,,=-4.o,^,=-4.5. (1 1) 

With increasing iron abundance, the da index increases, whereas 
the influence of silicon on 6a is only marginal. On the other hand, 
both the silicon and iron abundances influence the value of the 
6(h - y) index. This impUes that for a star with a given effec- 
tive temperature (and an assumed homogeneous surface distri- 
bution of elements), it is possible to infer the value of iron and 
silicon abundance directly from photometry. One has to keep in 
mind, however, that other elements may also influence the value 
of6ib-y) and 6a (cf., Kupka et al. 2004; Khan & Shulyak 2007). 
Consequently, such procedure is likely possible only for hot CP 
stars. Finally, for stars with a strong magnetic field, the Zeeman 
splitting of iron lines may also influence the a index (Khan & 
Shulyak 2006). 

As the star rotates, surface regions with difl'erent chemical 
composition appear on the visible disc and consequently 6a and 
6(b - y) also vary with rotational phase (see Figs. 11, 15). Here 
we plot the light curve calculated as described in Sect. 3.3 with 
filter parameters taken from Tab. 3. From this plot, we can con- 
clude that the amplitudes of 6a variations are very small, on the 
order of miUimagnitudes. 

6.5. Hp line profile variations 

The fact that we are able to predict the correct amplitudes of 
brightness variations in both the UV and visual spectral domains 
from just the abundance maps of iron and sihcon supports the 
conclusion that the effective temperature is constant on the sur- 
face of HR 7224. This finding weakens one of possible expla- 
nations of Balmer lines' profile variations as summarised by 
Lehmann et al. (2007). Since the effect of the Lorenz force on 
the line profiles is mainly on the fine wings (cf., Shulyak et al. 
2007), it cannot explain the line profile variabiUty observed in 
HR 7224. We tested if our model is also able to explain the ob- 
served Hj6 line variability. For this purpose we calculated the 
line profiles in individual phases using the code developed by 
Skalicky (2008). The code numerically integrates the emergent 
intensity IdO, Q) over the visible stellar surface Q. taking into ac- 
count the Doppler shift due to the stellar rotation. The emergent 
intensity from each surface point is calculated by interpolating 
the synthetic spectra I{A, 0, esi, spe) for appropriate surface abun- 
dance distribution of sihcon and iron (Lehmann et al. 2007). We 
assume fixed effective temperature and surface gravity. The pre- 
dicted phases of the maximum line depth (Fig. 16) roughly cor- 
respond to the observed ones (Lehmann et al. 2006, Fig. (13) 



therein), also the shape of the minimum at the phase <p « 0.2 
is roughly correct, however the observed amplitude of the min- 
imum at the phase (f> « 0.6 seems to be by a factor of about 2 
higher. We conclude that either another elements has to play the 
role in the Ufi line variability, or the observed difference between 
observation and theory is partially caused by somehow noisy ob- 
servational data. 



7. Conclusions 

We successfully simulated the light variabihty of HR 7224 di- 
rectly from the silicon and iron surface abundance maps derived 
by Lehmaim et al. (2007). There is very good agreement be- 
tween the observed and predicted light variability in the uvby 
bands of the Stromgren photometric system. We did not intro- 
duce any free parameter to improve the agreement between the 
theoretical and observed hght curves. 

The rotationally modulated light variability of HR 7224 is 
caused by the flux redistribution due to iron line transitions, sil- 
icon bound-free transitions, and by the inhomogeneous surface 
distribution of these elements. This picture is also supported by 
the agreement between the predicted behaviour of the UV flux 
distribution and that observed by the TDl satellite. 

We support the conclusion of Khan & Shulyak (2007) that 
numerous iron lines contribute significantly to the weU-known 
depression at 5200 A. Moreover, we show that iron is able to 
influence the peculiarity index a, whereas silicon's influence 
is marginal. With a suitable calibration, this could enable the 
derivation of abundances in hot CP stars directly from photome- 
try. 

We conclude that a promising explanation for the light vari- 
ations in CP stars is a flux redistribution through line and bound- 
free transitions combined with the inhomogeneous surface dis- 
tribution of various elements. 
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